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Abstract. I review particle heating by MHD turbulence in collision- 
less plasmas appropriate to advection-dominated accretion flows. These 
considerations suggest that the preferential turbulent heating of protons 
assumed by theoretical models is only achieved for relatively subthermal 
magnetic fields. 



1. Introduction 

A number of authors have argued that, at sub-Eddington accretion rates, the 
gravitational potential energy released by turbulent stresses in an accretion flow 
may be stored as thermal energy, rather than radiated as in thin accretion disks 
(Ichimaru 1977; Rees et al. 1982; Narayan k Yi 1994, 1995; Abramowicz et al. 
1995; see Narayan, Mahadevan, & Quataert 1998 and Kato, Fukue, & Mineshige 
1998 for reviews). In this case, the gas heats up to nearly virial temperatures 
and is assumed to adopt a two-temperature configuration, with the protons 
significantly hotter than the radiating electrons. 

There are two crucial microphysical issues relevant to these accretion mod- 
els, which are called either ion tori (Rees et al. 1982) or advection-dominated 
accretion flows (ADAFs; Narayan & Yi 1994). The first is the value of 5, the 
fraction of the turbulent energy in the plasma which heats the electrons (a frac- 
tion 1 — S heats the protons). The second is whether or not the dominant 
process for exchanging energy between electrons and protons is Coulomb colli- 
sions. ADAF models typically assume that <5 1 (preferential proton heating) 
and that Coulomb collisions represent the only thermal coupling. In this case, 
the plasma is unable to cool because it is nearly collisionless. The protons, 
which by assumption receive the turbulent energy, are (1) unable to radiate and 
(2) unable to transfer their thermal energy to the electrons, the more efficient 
radiators. 

The purpose of this review is to present some rather simple, but general, 
theoretical considerations regarding particle heating in ADAFs, which aim to 
address the value of 5. 
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2. Physical Picture 



A fundamental assumption of my analysis is that the particle distribution func- 
tions, while not necessarily strictly thermal, are not pathologically nonthermal; 
that is, I assume that the plasma does not reach an equilibrium configuration 
with, e.g., highly anisotropic or bump- in-tail distribution functions. 

For the simple considerations presented here, only two properties of the 
plasma are important. The first is the proton to electron temperature ratio, 
T p /T e , taken to be 3> 1. Note that it is consistent to assume T p » T e while si- 
multaneously investigating particle heating. This is because the two temperature 
nature of the plasma is in large part due to the efficiency with which relativistic 
electrons cool; it is rather insensitive to 5. The second property of importance 
is the magnetic field strength in the plasma, parameterized by f3 = P g as/-fmag- 
The quantity f3 must be > 1 for the gravitationally confined plasmas of interest. 

Like thin disks, ADAFs are turbulent magnetized plasmas. On scales much 
larger than the Larmor radius of the thermal protons (= p), we can treat the 
turbulence using MHD, superposition of the magnetosonic and Alfven 

waves. Particles are heated when the turbulence cascades to sufficiently small 
scales that it is dissipated. 

3. Collisionless Damping 

Let a wave in the plasma have wave vector k and frequency u and define _L and 
|| to be directions perpendicular and parallel to the local magnetic field. 

In collisionless plasmas, the wave dissipation mechanisms of interest are 
wave-particle resonances (molecular viscosity, thermal conductivity, electrical 
resistivity, etc. are entirely unimportant). In the next section, I will argue that 
the turbulence always has frequencies well below the proton cyclotron frequency. 
In this case, resonance occurs when the wave's phase speed along the field line, 
v = cj/fe||, equals v», the velocity of a particle along the field line. A necessary 
(but not sufficient) condition for strong damping is that v be comparable to 
the thermal speed of the particles, so that there are a large number of resonant 
particles. 

This resonance actually corresponds to two physically distinct wave-particle 
interactions. In Landau damping (LD), particle acceleration is due to the wave's 
longitudinal electric field perturbation (i.e., the usual electrostatic force, Eh). 
For a wave damped solely by LD, 6 ~ 1 since the electrons are preferentially 
heated by a factor of w (m e T p 3 /rripT^) 1 / 2 , which is > 1 for T p ~ 10 2 - 10 3 T e . In 
transit-time damping (TTD), the magnetic analogue of LD (and the collisionless 
analogue of Fermi acceleration; Achterberg 1981), the interaction is between the 
particle's effective magnetic moment (fi = mv\/2B) and the wave's longitudi- 
nal magnetic field perturbation, Bu (Stix 1992). For a wave (with k^p < 1) 
damped solely by TTD, the protons are preferentially heated by a factor of 
w {m p T, p /m e T e ) 1 ^ 2 3> 1; the contribution to 5 is therefore always <C 1. This is 
because in plasmas with T p 2> T e , the protons have the larger magnetic moment 
and so couple better to a wave's magnetic field perturbation. 



2 



4. MHD Turbulence 



In a plasma with > 1, one can (roughly speaking) decompose MHD turbulence 
into three modes, the Alfven wave and the fast and slow magnetosonic modes. 
The fast mode has the dispersion relation uj « ^ p |k| (v p is the proton thermal 
speed) and is essentially a sound wave. The slow mode has oo « (va ~ 

VpP" 1 / 2 is the Alfven speed) and, for kj_ ^ 0, has a parallel magnetic field 
perturbation (-By). The Alfven wave also has to ~ but is incompressible, 

with £j| = 5|| = in the limit that the wavelength is long compared to the 
Larmor radius of thermal protons (|k|p 1). 

4.1. Sound Waves (Fast Modes) 

Sound waves in an unmagnetized collisionless plasma are strongly damped by 
LD. For T p » T e , the dissipation time is formally <C the mode period, implying 
that the modes are non-propagating. These modes are unlikely to be efficiently 
excited since, quite generally, strongly damped oscillators dissipate less energy 
than weakly damped oscillators (consider, e.g., the analogous problem of two 
resistors in parallel; the resistor with the smallest resistance dissipates the most 
power). The same holds true for sound waves in a j3 > 1 plasma, provided 
/cj_ ~ fey. For k± 3> fey, i.e., nearly perpendicular sound waves, the parallel 
phase speed is v « v p k±/k\\ 3> v p , v e (since v e ~ v p by T p 3> T e ). Since the wave 
is highly suprathermal, there are very few particles for the wave to resonate with 
and it is consequently undamped. These waves likely dissipate by steepening 
and forming weak shock waves. 

In what follows I will neglect quasi-perpendicular sound waves, assuming 
that they are energetically unimportant. The reason is that I believe that the 
turbulence in ADAFs is intrinsically magnetic in origin, being initiated by some 
combination of Balbus-Hawley and convective instabilities. The turbulent speed 
(vt) is therefore ~ v a < v p . Moreover, I will argue that va is likely to be <C v p , 
in which case, if v t ~ va, excitation of sound waves will be strongly suppressed 
(subsonic instabilities do not easily excite sound waves). 

4.2. Slow Modes 

In a > 1 collisionless plasma, slow modes are strongly damped and non- 
propagating for fci / (Foote & Kulsrud 1979). The reason is that, even in 
the limit of wavelengths long compared to the Larmor radius of thermal protons 
(|k|p <C 1), slow modes have a parallel magnetic field perturbation (Bu ^ 0) and 
are damped by TTD. They £1X6 . QjS with sound waves, unlikely to be excited. For 
/c_i_ — ► 0, these considerations fail, since B\\ — > 0. Quasi-parallel slow modes are, 
however, likely to cascade in k±, develop By ^ 0, and damp by TTD (heating 
the protons). 

4.3. Alfven Waves 

The Alfven wave has a parallel phase speed \v\ = va', in plasmas appropriate 
to ADAFs, va is comparable to the electron and proton thermal speeds and so 
there are a large number of particles available to resonate with the wave. In 
the MHD limit, however, the Alfven wave is undamped by linear collisionless 
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effects. This is because En = and B» = 0, i.e., there is no electric/magnetic 
field perturbation associated with the wave which can accelerate particles. 

This lack of dissipation on large scales implies that understanding heating of 
particles by Alfvenic turbulence requires understanding how the energy cascades 
to small scales, where dissipation finally becomes important. 

Alfvenic Turbulence The nature of Alfvenic turbulence in a collisionless, (3 > 1, 
plasma is not fully understood. A plausible scenario is due to Goldreich and 
Sridhar (1995; GS). They argue that Alfvenic turbulence naturally evolves into 
a critically balanced state in which the timescale for nonlinear effects to transfer 
energy from a wavevector ~ k to a wavevector ~ 2k (= the cascade time, T c ) is 
comparable to the linear wave period at that scale, T = 2ituj~ 1 ; this determines 
how rapid the dissipation must be to halt the cascade. It also implies that the 
cascade is highly anisotropic, with the energy cascading primarily perpendicular 
to the local magnetic field; the parallel and perpendicular sizes of a wave at any 

scale are correlated, with fen ~ fc^ 3 L -1 / 3 <C k±, where L is the outer scale of 
the turbulence. 

Quasi-perpendicular Alfvenic turbulence has the important property that 
its frequencies are always much less than the proton cyclotron frequency. This is 
because the cascade occurs primarily in k± while the Alfven frequency is oc fey . 
This justifies my claim in §3 that low frequency wave-particle resonances are of 
primary importance. 

Alfven Wave Damping The character of an Alfven wave changes when k±p ~ 1; 
at this point the wave develops a non-zero parallel magnetic field perturbation 
(in contrast to long wavelength perturbations, for which B» = 0). The waves 
can therefore be dissipated by TTD (Quataert 1998, Gruzinov 1998, Quataert 
& Gruzinov 1998; QG). 

Figure 1 shows kinetic theory calculations of the dissipation of nearly per- 
pendicular Alfven waves. For k±p < 1 (where the waves can still be called 
Alfven waves), the dissipation is independent of T p /T e , but is a strong function 
of (3. The electron temperature is unimportant since the electrons do not partic- 
ipate significantly in the damping. The dissipation is sensitive to (3 ~ (v p /va) 2 , 
since this determines the number of resonant protons and the protons' magnetic 
moment (which is the effective coupling coefficient between the protons and the 
wave) . 

For large /3 > 1, Alfven waves are strongly damped for kj_p ~ 1. Most of 
the turbulent energy heats the protons, since they are responsible for TTD (§3). 
For (3 = 1, however, the maximal dissipation rate for an Alfven wave in the GS 
cascade (obtained at k±p = 1) is 7T 0.1 (see Figure 1). Since the timescale 
for energy to cascade through the inertial range is T c ps T, this suggests that, for 
plasmas with equipartition magnetic fields {(3 ~ 1), very little of the turbulent 
energy is dissipated on scales comparable to or greater than p. 

Alfven waves only exist for k±p < 1. For kj_p > 1, the same mode is 
called the whistler. As emphasized by Gruzinov (1998; see also QG), whistlers 
are unlikely to heat the protons. For k±p > 1, but k± 3> fey, whistlers have 
uj <C Q p . Thus, in a mode period, a particle undergoes many Larmor orbits. 
Since the mode's perpendicular wavelength is smaller than the proton Larmor 
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Figure 1. Dimensionless dissipation rate of an Alfven wave (kj_p < 1) 
or whistler wave (k±p > 1) as a function of perpendicular wavenumber. 

radius (kj_p > 1), the protons (but not the electrons) sample a rapidly varying 
electro-magnetic field in the course of a Larmor orbit. As a result, they are 
"frozen out" and become dynamically unimportant; the protons effectively no 
longer "see" the wave. They therefore cannot contribute to damping the whistler 
energy, which cascades to smaller length scales until it is damped by the electrons 
(this is the origin of the strong damping at k±p 3> 1 in Figure 1). 

The Physics of 5 The above arguments suggest that the crucial physics de- 
termining the electron heating rate in ADAFs is what fraction of the turbulent 
energy is dissipated as Alfven waves (heating the protons) and what fraction 
cascades to small scales, becoming whistlers and heating the electrons. This 
depends on (3, but not on T p /T e . Unfortunately, estimates of the fraction of the 
turbulent energy that becomes whistlers are highly uncertain, as they are expo- 
nentially sensitive to the details of turbulence on the scale of the proton Larmor 
radius, where both nonlinear and kinetic effects are important. For example, we 
do not accurately know T c , the nonlinear timescale of the turbulence. Nor is the 
calculation of the damping of the turbulence that secure (because it is based on 
a theory of linear Alfven waves damped by thermal particles; both linearity and 
thermal protons are crude approximations). 

Parameterizing the uncertainty in the details of the turbulence, and using 
numerical simulations of MHD turbulence as a guide, QG estimated that the 
critical (3 above which proton heating dominates (say 5 < 0.1) is ~ 10, with an 
uncertainty of ~ 10! 

5. Summary and Future Work 

The physics of particle heating in a collisionless plasma suggests that, in ADAFs 
and ion-tori, proton heating may be favored for relatively subthermal magnetic 
fields, but not for strictly equipartition ones {{3=1). 
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In this paper, I have focused primarily on heating of particles by small scale, 
incompressible, Alfvenic turbulence. For (5 S> 1, it should be quite accurate to 
neglect sound wave and slow mode excitation. As emphasized by Blackman 
(1998), however, one expects the importance of compressibility to be an in- 
creasing function of decreasing (5. For example, at (3 = 1, slow mode excitation 
cannot be neglected (these modes are damped by TTD on relatively large scales, 
heating the protons; e.g., Blackman 1998). Compressibility (and hence dissipa- 
tion) of large scale Alfven waves is also potentially important. My belief is that, 
while these effects are undoubtedly present, it is unlikely that the majority of 
the turbulent energy is dissipated on large scales; thus, to order of magnitude, 
the analysis presented in this review is appropriate. More quantitative estimates 
are, however, clearly needed. 

The unknown importance of magnetic reconnection, and its (presumed) 
electron heating, remains a significant uncertainty for ADAF models (Bisnovatyi- 
Kogan & Lovelace 1997; see also QG). 
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